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We propose a model to explain the reduced bimolecular recombination rate found in state-of-the-art bulk 
heteroj unction solar cells. When compared to the Langevin recombination, the experimentally observed rate 
is one to four orders of magnitude lower, but gets closer to the Langevin case for low temperatures. Our 
model considers the organic solar cell as device with carrier concentration gradients, which form due to the 
electrode/blend/electrode device configuration. The resulting electron concentration under working conditions 
of a solar cell is higher at the cathode than at the anode, and vice versa for holes. Therefore, the spatially 
dependent bimolecular recombination rate, proportional to the local product of electron and hole concentration, 
is much lower as compared to the calculation of the recombination rate based on the extracted and thus averaged 
charge carrier concentrations. We consider also the temperature dependence of the recombination rate, which 
can for the first time be described with our model. 

PACS numbers: 7 1.23. An, 72.20Jv, 72.80.Le, 73.50.Pz, 73.63.Bd 
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I. INTRODUCTION 



Organic bulk heteroj unction solar cells have shown an in- 
creasing performance in the recent years, and also scien- 
tific progress concerning the fundamental understanding has 
been made. 12 However, the dominant charge carrier loss 
mechanism determining the photocurrent is still under dis- 
cussion. The relevant processes during which the losses 
can occur are geminate recombination during polaron pair 
dissociation 34 , nongeminate recombination during transport 
of the already separated polarons 56 , and charge extraction 
from the device.™ A detailed analysis considering the inter- 
play of these mechanisms has still to be done, as already the 
separate processes are not completely described yet. Concern- 
ing nongeminate recombination, classically the bimolecular 
Langevin formalism 9 10 has been used for low mobility mate- 
rials. However, already in 1997, reports on a reduced rate as 
compared to Langevin 's derivation were discussed for polaron 
recombination in conjugated polymers, 11 and recently a sim- 
ilar reduction was found for polymer-fullerene solar cells P21 
Last year, we presented investigations of the polaron recom- 
bination in poly(3-hexyl thiophene):[6,6]-phenyl-C61 butyric 
acid methyl ester devices. We found a reduced recombination 
rate as compared to classic Langevin recombination, with a bi- 
molecular decay in pristine samples, and a third order recom- 
bination in annealed samples. 6 The low recombination rate 
as well as the third order decay have been observed by other 
researchers as well, 5 13 the origin of both effects remaining 
unresolved. For the third order recombination, the scenario 
of a carrier concentration or time dependent bimolecular re- 
combination should be considered, and is probably related to 
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delayed recombination due to trapping in the tail of the den- 
sity of stat es!^ Concerning the reduced recombination rate, 
two mode kr^^ have been proposed in literature trying to ex- 
plain the reduction mechanism. However, as we pointed out 
recently, 6 both fail to predict the correct temperature depe- 
nence. In this paper, we will present a simple model predict- 
ing the low bimolecular recombination rate as compared to the 
Langevin theory, as well as its temperature dependence. 



II. MODEL 

A. Existing models for the reduced Langevin recombination 

Before introducing our model explaining the reduced 
Langevin recombination, let us briefly present the basic idea 
behind the previously published models. Four different re- 
combination models are shown in Fig. [TJ (a) the classic 
Langevin recombination, (b) the minimum mobility model by 
Koster et alP3 (c) the model by Adriaenssens and Arkhipov, 11 
and (d) our model. 

In the classic Langevin recombination (Fig. [TJa)), the 
derivation of which is nicely shown in the book by Pope and 
Swenberg, 10 assumes that the rate limiting factor for recombi- 
nation is the finding of the respective recombination partners 
(1), and not the actual recombination rate (2). Neglecting pro- 
cess (2) as it is faster than (1), the finding of electron and 
hole depends on the sum of their diffusivities or — considering 
the Einstein relation — their mobilities. Thus, the Langevin re- 
combination rate is 

R = y(np-nj), (1) 
where n and p are electron and hole concentrations, respec- 
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FIG. 1 : (Color Online) Recombination mechanisms in low mobility 
materials, all being based on the Langevin recombination shown in 
(a), (a) to (c) consider local regions within a device, whereas (d) 
corresponds to the whole device of thickness L. n(x) and p(x) are 
the position dependent electron and hole concentration, respectively. 
The details are described in the text. 



tively, and ni is the intrinsic carrier concentration. Here, 

-(jUe+JLlh) (2) 



£ r £0 



is the Langevin recombination prefactor, where q is the el- 
ementary charge, 8 r 8o the effective dielectric constant of the 
ambipolar semiconductor, ju e and jiih the electron and hole mo- 
bilities. 

The model introduced by Koster et al!^ is an extension 
of the Langevin model. It considers an inherent property of 
bulk heteroj unction solar cells: the phase separation of donor 
and acceptor materials (Fig. [TJb)). Under the reasonable as- 
sumption that holes are exclusively transported in the donor 
polymer phase, and electrons through the fullerene acceptor, 
a bimolecular recombination can only take place at the het- 
eroj unction. Therefore, if the slower charge carrier does not 
reach the interface, no recombination takes place. In order 
to consider this behaviour, Koster et al. let the recombination 
prefactor be governed by the minimum mobility, 



£ r £0 



mm(ju e ,ju h ). 



(3) 



The Arkhipov model 11 proposes that potential fluctuations 
in an ambipolar material are responsible for the recombination 
rate reduction (Fig. [IJc)). As the band gap remains constant, 
electrons and holes accumulate at the potential minima of the 
corresponding bands, therefore being spatially separated. In 
order to recombine, a potential barrier proportional to the en- 
ergy difference between the minimum and maximum of the 
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FIG. 2: (Color Online) Langevin recombination reduction factor £ in 
dependence on temperature. Shown are the results for annealed de- 
vices of ours (diamonds) P and of Juska et alP^ (circles). The models 
by Koster et alPl (dashed) and Arkhipov et alPl (dotted) are also in- 
cluded; they show a markedly different temperature dependence as 
compared to the experimental data. 



band fluctuations has to be overcome. In some respect, this 
model is similar to the Koster model, in as far as it also ac- 
counts for a spatial separation of the recombination partners. 
The origin of such a spatial separation could also be due to 
the above mentioned donor-acceptor phase separation. The 
recombination prefactor in the framework of the Arkhipov 
model is changed to 
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exp 



-«»), 



(4) 



where AE is the activation energy, and kT the thermal energy. 
If the Arkhipov model could be applied to bulk heteroj unc- 
tion solar cells, one would expect AE to be proportional to the 
energy difference between either the polymer and fullerene 
lowest unoccupied molecular orbitals, or between the poly- 
mer and fullerene highest occupied molecular orbitals. 

In our recent publication on bimolecular recombination, 6 
we simultaneously determined the time dependent carrier con- 
centration and charge carrier mobility in a photoinduced CE- 
LIV (charge extraction by linearly increasing voltage) experi- 
ment on pristine and annealed poly(3-hexyl thiophene):[6,6]- 
phenyl-C61 butyric acid methyl ester solar cells. A positive 
aspect of this experimental technique is that is able to deter- 
mine the reduction factor 



R 



R, 



(5) 



experiment 



which is the fraction of the Langevin rate R (Eqn. ([T])) over the 
experimentally determined recombination rate ^experiment- A 
similar experiment had been done previously by Juska et alPS 
From the experimental results shown in Fig. [2] £ shows a neg- 
ative temperature coefficient in annealed devices, increasing 
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with decreasing temperature, thus minimizing the difference 
to the classic Langevin rate. In constrast to the experimen- 
tal findings, the Arkhipov model as well as the Koster model 
predict a positive temperature coefficient of thus increasing 
with temperature, going asymptotically closer to the classic 
Langevin rate = 1). As a side note, in order to calculate 
the temperature dependence for the Koster model, tempera- 
ture dependent mobilities have to be used. A suitable model 
is the gaussian disorder modelj^ which implies an exponen- 
tially decreasing mobility with falling temperature, depending 
mostly on the energetic width of the gaussian density of states, 
c, also called disorder parameter. Attributing different values 
of a to the electron and hole transporting phases in a bulk het- 
erojunction, which is in line with Koster's assumptions, the 
temperature dependence of Eqn. ^ shows the behaviour as 
described above. Only in the case of having the same disorder 
parameter for electrons and holes does the calculated recom- 
bination prefactor Jk become temperature independent, but it 
never can attain the experimentally found temperature coeffi- 
cient. Thus, neither the Juska^ nor the Koster 15 model can 
predict the temperature dependence of the recombination re- 
duction factor £ correctly. 



B. Carrier concentration gradient model for the reduced 
Langevin recombination 

Our model considers the discrepancy between the experi- 
mental determination of the carrier concentration by charge 
extraction techniques, which gives only average values, and 
the locally varying carrier concentration gradients found in the 
devices under working conditions. This difference is of par- 
ticular importance for bulk heteroj unction devices, which con- 
sist of an ambipolar semiconductor layer — the donor-acceptor 
blend — sandwiched between anode and cathode. 

The experimental carrier concentration yields the complete 
density of charge carriers ^extracted found in the device under 
test conditions. When calculating the Langevin recombina- 
tion rate, Eqn. ([I]), usually the assumption ^extracted = n = p is 
made, so that R = Y^extracted should fit the experimental data. 
However, as described above, the additional reduction factor 
£ had to be introduced in order to yield a good description of 
the experimental carrier concentration decay with time, with 
^experiment = ^extracted bein S the experimentally determined 
recombination rate. 

The problem with the above mentioned assumption of 
^extracted = n = p stems from the implicit consequence n(x) = 
p(x), where x is the distance from anode to cathode of the 
device. Even if electrons and holes could be extracted sepa- 
rately, only spatial averages n(x) and p(x) were experimen- 
tally accessible. However, considering the charge carrier dis- 
tribution in a bulk heteroj unction device, these are not valid 
assumptions under most measurement conditions. A sketch of 
a typical electron and hole carrier concentration in a bulk het- 
erojunction solar cell is shown in Fig.[TJd). The strong carrier 
concentration gradients are indeed typical for an ambipolar 
device with asymmetric contacts such as a bulk heteroj unc- 
tion solar cells. In principle, these gradients occur in the dark 



case and under illumination; in the latter case, the carrier con- 
centration gradient is somewhat lower due to the photogener- 
ation of electron-hole pairs, but nevertheless very relevant to 
the topic under discussion. The limiting factor for the recom- 
bination is still the finding of electron and hole (1), which is 
proportional to the sum of the mobilities, as described above. 
However, now the different electron and hole concentration 
gradients have to be considered. The indium tin oxide (ITO) 
electrode, the anode, is a good hole injection contact into con- 
jugated polymers such as P3HT. Consequently, the hole con- 
centration of the whole device finds its maximum at this spa- 
tial position for voltages below the built-in voltage. At the 
same time, hole concentration at the cathode is much lower 
for voltages below the flat band case, in darkness and under 
illumination. The concentration gradients are lowered due to 
illumination, as the generation of electron-hole pairs through- 
out the bulk changes the carrier concentration mostly where it 
was low without light. That means that the relative increase 
of the hole concentration is highest in the vicinity of the cath- 
ode, despite the extraction path of the photogenerated holes 
being via the anode. Thus, the considerations concerning the 
carrier concentration gradients apply to both, dark and illumi- 
nated devices. In order to calculat e a re com binat ion rate based 
on average carrier concentrations n(x) and p(x) — considering 
these conditions — introduces a large systematic error. 

In order to better illustrate the differences arising when 
comparing recombination rates calculated from either aver- 
age carrier concentrations or actual gradients, we devised a 
simple model. This model will allow us to get a better im- 
pression of the recombination reduction factor and thus the 
origin of the seemingly reduced Langevin recombination rates 
in ambipolar organic devices, in particular bulk heteroj unction 
solar cells. 

According to our statement, £ can be defined as 

r zfo n ( x )p( x ) dx ,~ 

^ = ~T\ ~T\ ' (6) 

n[x) - pyx) 

where the denominator corresponds to spatial averages, for in- 
stance when using carrier concentrations from charge extrac- 
tion experiments. The numerator instead correctly accounts 
for the carrier concentration gradients found in the device un- 
der a certain applied voltage and given light intensity. 

For simplicity, we define mirror- symmetric electron and 
hole carrier concentration gradients by 

n(x) = n n exp(-a^) (7) 
p(x) = p p Qxp (^-a^j-^j . (8) 

Here, a = \n{n n /n p ) = hx(p p /p n ), where n n (p p ) is the elec- 
tron (hole) concentration at the electron (hole) injecting con- 
tact, and n p (p n ) is the electron (hole) concentration at the 
anode (cathode), x is the distance from anode (0) to cathode 
(L). The resulting distance dependent carrier concentrations 
are shown in the inset of Fig. [3] 

Now, we can continue the calculation started in Eqn. ^ 
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FIG. 3: (Color Online) The simple model for the Langevin reduc- 
tion factor given by Eqn. |To| . It depends on the parameter a, 
which represents the carrier concentration gradients. The inset shows 
the carrier concentration (on a logarithmic scale) vs. distance after 
Eqns. (|7J and ([8]). p p (n p )is the hole (electron) concentration at the 
hole injection electrode, the anode. Similarly, n n (p n )is the electron 
(hole) concentration at the electron injection electrode, the cathode. 

using Eqns. ([7]) and ([5]), 

^ = zfo n ( x )p( x ) dx (9) 
zfo n ( x ) dx 'zfoP( x ) dx 
= a 2 e- a (l-e- a y 2 (10) 

The calculated Langevin recombination reduction factor is 
shown in Fig. [3] The steeper the carrier concentration gra- 
dients, the larger the discrepancy to their respective spatial 
averages, the smaller Two important consequences arise: 
first, the recombination reduction factor £ depends on all pa- 
rameters changing the carrier concentration gradients, such as 
applied voltage and charge carrier mobility. Second, £ does 
not directly depend on the actual recombination mechanism: 
no matter if the dominant recombination is monomolecular or 
bimolecular, or if the solar cell is illumated or not, what counts 
is the resulting steady state carrier concentration. 

In macroscopic device simulators considering at least one 
space dimension, carrier concentration gradients are already 
accounted for. Therefore, £ does not need to be explicitly con- 
sidered. On the contrary, as the simple model presented above 
uses very rudimentary functions to describe the carrier con- 
centration gradients, not considering injection barriers etc., 
we will apply a macroscopic simulation program in order to 
better understand the apparently reduced recombination under 
typical measurement conditions. 

III. MACROSCOPIC SIMULATION 

The macroscopic simulation program implemented by us 
solves the differential equation system of the Poisson, conti- 



parameter value description 

E Gap 1 . 1 eV effective bandgap 1 ^ 

<&p 0.1 eV injection barriers 

/u n , n p 1 • 10~ 8 m 2 V~ 1 s- 1 mobilities 

d 100 nm active layer thickness 

G 6.0- 10 27 m _3 s _1 generation rate 

T 300 K temperature 

N e ff 1 .0 • 10 26 m -3 effective density of states 

8 r 3.4 relative static permittivit) 

TABLE I: Parameters used in the macroscopic simulation. 

nuity and drift-diffusion equations by an iterative approach as 
described in Ref. 8 . Additionally, we consider injection barri- 
ers at both electrodes as well as a finite surface recombination. 
We use the field independent surface recombination velocity 
S(0) of the well known Scott-Malliaras model 22 considering 
mirror charge effects at surfaces for both electrodes. The sur- 
face recombination current is defined as 

jRec = qS(0)(n-n th ) (11) 

with 

5(0) = 16tt88 (kT) 2 /u/q 3 . (12) 

Here, n is the electron concentration at the surface, and n tn 
is the thermally activated carrier concentration. The surface 
recombination current is defined for both carrier types at each 
electrode. 

In order to clearly and unambiguously show the effect of 
internal charge carrier distribution and its impact on the re- 
combination rate, we deactivated the field dependent polaron 
pair dissociation. Consequently, the net generation rate U (x) 
is simply a function of the generation rate G and the classical 
Langevin recombination as defined in Eqn. ([T]), 

U(x) = G— (/nn {n(x)p(x)-n}) . (13) 

The parameters assembled in Tab. [I] were used for all simu- 
lations, unless explicitly mentioned. Where temperature de- 
pendent calculations were performed, we varied the mobility 
according to the gaussian disorder model 17 with a disorder 
parameter of a = 75 meV and a prefactor chosen to achieve a 
mobility of 10 -8 m 2 /Vs for electrons and holes at 300K. 

We point out that macroscopic simulations are very use- 
ful to study organic devices such organic bulk heteroj unction 
solar cells, despite the assumption of an effective medium. 
In the latter, the hole conducting properties are derived form 
the donor material, whereas the electron conduction properties 
come from the acceptor material. For donor-acceptor blends 
with a very fine-grained phase separation, the assumption of 
an effective medium is very good. For coarser phase separa- 
tions, the situation becomes more difficult, as band bending 
between the two phases cannot be described with the effective 
medium. Nevertheless, as parameters derived from micro- 
scopic Monte Carlo and Master equation simulations 23 2425 



5 




20 40 60 80 100 
Distance from Anode [nm] 



FIG. 4: (Color Online) Simulated steady state electron and hole con- 
centrations in bulk heteroj unction solar cells under short circuit (top 
graph), open circuit (middle), and at the built-in potential (bottom) at 
an illumination of 1 sun. In each graph, the concentration profiles for 
the balanced electron-hole mobilities \x = 10 -5 (short dashed line), 
10 -8 (solid line), and 10 -11 m 2 /Vs (long dashed line) are shown. 
Holes have the highest concentration at the anode (left), electrons at 
the cathode (right). 



and analytic theory 26 27 can be used to describe the properties 
of the donor-acceptor blend, the use of macrosopic simula- 
tions offers a very good insight into the impact of microscopic 
charge transport and recombination properties on the macro- 
scopic device parameters such as current-voltage character- 
istics. Indeed, macroscopic simulations complement the mi- 
croscopic point of view very well, in particular as also asym- 
metries due to different work functions for electron and hole 
injection, and their influence on the device properties can be 
studied. The usefulness of this approach has been reported 
previously ! 8 * 28 ! 29 ! and is in the focus of the present work as 
well. 



IV. RESULTS AND DISCUSSION 

Fig. [4] shows the electron and hole concentrations under 
short circuit, open circuit, and the built-in potential. The latter 
is the voltage at which photo-CELIV measurements are usu- 
ally performed. In comparison to our simple model (Fig. [3]), 
the concentration gradients have a more complicated shape 



due to injection and extraction as well as the interplay of gen- 
eration and recombination in steady state. Nevertheless, it is 
clear that the effect remains the same. For each case shown, 
£ as defined in Eqn. ^ is much smaller than unity. Conse- 
quently, charge extraction experiments such as photo-CELIV 
or transient photocurrents will yield the average carrier con- 
centrations, which — if used to calculate the bimolecular re- 
combination rates — will yield overestimated values. Looking 
in more detail, some important features of £ are seen. The 
carrier concentration gradient is strongest at short circuit, and 
very low at the built in potential. In the former case, charge 
extraction is most favourable, whereas under flatband condi- 
tions, the charges tend to stay within the device due to the 
lack of a driving force. Similarly, a high mobility will tend to 
create steeper gradients. 

A more general feature of the bimolecular recombination 
in an ambipolar device is implicitly shown in Fig. |4j due to 
the opposite electron and hole concentration profiles as well 
as the influence of the electrodes, the local product of electron 
and hole density is lowest where the deviation of the both con- 
centrations to the average is largest — and correspondingly, £ 
is smallest. This means that the polaron recombination is gen- 
erally weakest at the contacts. 

The discrepancy between average and local product of elec- 
tron and hole concentrations in dependence on the charge car- 
rier mobility is shown in Fig. [5jtop), the resulting recombi- 
nation reduction factor £ in Fig. [5jmiddle). Both, mobili- 
ties and injection barriers, were chosen to be symmetric for 
electrons and holes, but qualitatively the results hold true for 
asymmetric conditions as well. A high mobility corresponds 
to an efficient charge extraction, leading to steeper electron 
and hole concentration gradients, and thus a lower Simi- 
larly, the lower the injection barrier, the weaker the concen- 
tration gradients, the closer is £ to unity. The corresponding 
bimolecular recombination rates are shown in Fig. [^bottom). 
Also included is R/C, (Eqn. ([5])), the recombination rate as de- 
rived when only considering average carrier concentrations. 
It equals yn 2 , thus implying a severe overestimate of the loss 
rate. 

Using our macroscopic simulation, we find two other influ- 
ences as well. Raising the external voltage from zero to the 
built-in voltage, £ will gradually converge to unity, as the car- 
rier concentration gradients become level when coming closer 
to flat band conditions. The effect of photon absorption is sim- 
ilar: electron-hole pair generation throughout the extent of the 
device leads to weaker carrier concentration gradients; conse- 
quently, the recombination reduction factor £ approaches one 
for high illumination densities (not shown). 

Photo-CELIV is an experimental technique which is able to 
determine the correct recombination rate present in the device, 
despite the fact that it only considers averaged recombination 
rates. This is possible due to the direct fitting of the experi- 
mental time dependent carrier concentration data to the charge 
carrier continuity equation, 

£ = (14) 
with the charge carrier generation rate G(t), which is zero for 
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FIG. 5: (Color Online) (Top) Squared carrier concentrations, by 
global (up) or local multiplication (rip) of electron and hole concen- 
tration, respectively, in dependence on the charge carrier mobility. 
The data shown was calculated for flat band conditions, as they are 
typically used in photo-CELIV measurements. (Middle) The recom- 
bination reduction factor in dependence on the charge carrier mo- 
bility, for different injection barriers (& p at the anode, <3? n at 
cathode) of eV (short dashed line), 0.1 eV (solid line), and 0.2 eV 
(long dashed line). £ is lower for efficient extraction at high mo- 
bilities, and for low injection barriers. (Bottom) The simulated bi- 
molecular recombination rate vs. the charge carrier mobility, R(jn), 
after Eqn. |TJ, as a result of the carrier concentration gradients. Also 
shown is the recombination rate which would be determined if 
the carrier concentration gradients were neglected. 



time t > due to using a nanosecond laser pulse, the recombi- 
nation reduction factor and the Langevin recombination rate 
R (Eqn. ([T]). Monomolecular contributions are neglected. As 
photo-CELIV yields n and jn simultaneously, R is completely 
known, and £ can be determined. 

For a comparison of our macroscopic simulation to experi- 
mental photo-CELIV data, see Fig. [6] Both data and simula- 
tion are for flat band conditions, i.e., the built-in potential. The 
shape of the experimental data of RefsP^is very well repro- 
duced by the simulations. As pointed out above, the models by 
Koster et all^and Adriaenssens and Arkhipov 11 were not able 
to describe this temperature dependence. Thus, our model is 
the first one to describe the experimentally found temperature 
dependence of the reduced Langevin recombination qualita- 
tively. Additionally, the voltage dependent £ as determined 
by Juska et alP^ corresponds to our simulations (not shown). 
We point out that the absolute magnitudes of simulated and 
experimentally determined reduction factors differs by a fac- 
tor around 0.1 to 0.005, the discrepancy being independent of 
electric field and temperature. Thus, the recombination reduc- 
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FIG. 6: (Color Online) Comparison of the temperature dependent 
recombination reduction factor £ for the experimental photo-CELIV 
data (as already shown in Fig. [2) with our macroscopic simulation. 
The difference of the absolute values of £ in simulation and experi- 
ment, a factor of 1/200 for the data of Juska et al. (circles) compares 
do the calculated values with an injection barrier of 0.1 eV (black 
solid line). The data of Deibel et al. (diamonds) has similar shape, 
with an additional, static reduction factor of about 1/20 is needed to 
match the simulation. The details are described in the text. 



tion factor £ is composed of two contributions, 

£(!T, F, G) = Cgradient(?\ F, G) • £ sta ti 



(15) 



The first term is the temperature T, electric field F, and charge 
carrier generation rate G dependent prefactor, which is due to 
the carrier concentration profiles in the device, as described by 
our model. The second contribution to £ is constant, and not 
considered in our simulation. This static contribution Static 
can be due one or more of the following factors: (a) a geomet- 
rical factor due to the donor-acceptor phase separation, the 
charges being confined to their respective phase, (b) the donor 
resp. acceptor material of the polymer-fullerene blend can 
have different dielectric constants, 30 or (c) deviations from the 
Langevin recombination factor due to energetic disorder, size 
of the donor-acceptor domains, and mismatch between the 
electron and hole mobility!^ 

Szmytkowski 30 calculates a temperature independent re- 
combination reduction factor, 



(16) 



with 8d being the relative permittivity of the donor, and 8 a 
the corresponding value for the acceptor material. For P3HT 
and PCBM, the relative permittivities are approximately 3.4 
and 4.0, respectively, thus ^ £ = 0.08, and even smaller if the 
respective permittivities are closer to one another. However, to 
our knowledge, this explanation has not been experimentally 
verified as of yet. 

Groves et al. 25 perform Monte Carlo simulations of elec- 
trons and holes in a blend system with hopping transport, and 
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study deviations from the Langevin theory of bimolecular re- 
combination inin view of bulk devices and field effect transis- 
tors. For bulk heteroj unctions, they point out that the effect of 
energetic disorder, domain sizes and electron-hole mobility 
mismatch leads to ^static of only between 0.1 and 1. There- 
fore, the authors suggest to consider the influence of deep car- 
rier trapping for explaining smaller recombination reduction 
factors. 

Thus, our model considering the carrier concentration gra- 
dients of electrons and holes in an ambipolar organic device 
can explain the behaviour of the reduced Langevin recombi- 
nation in terms of temperature and electric field dependence, 
an additional static contribution is needed to match the exper- 
imental recombination rates. 



V. CONCLUSIONS 

In conclusion, we have presented a simple model describ- 
ing the reduced Langevin recombination in organic solar cells. 
The origin of the reduction factor is based on two contri- 
butions: first, the discrepancy between average electron and 
hole concentrations considered in charge extraction experi- 
ment and the usually steep carrier concentration profiles in 



organic semiconductors. Second, a static factor related to the 
phase separation of donor and acceptor, their energetic dis- 
order and relative permittivities, the origin of which is still 
under discussion. Concerning the first contribution, the spa- 
tially dependent bimolecular recombination rate, proportional 
to the local product of electron and hole density, is much lower 
as compared to the rate based on average charge carrier con- 
centrations. The latter leads to an overestimation of the re- 
combination rates. Our model for the first time correctly de- 
scribes the qualitative temperature dependence of the reduc- 
tion factor found experimentally by the photo-CELIV method 
in P3HT:PCBM solar cells. It is also applicable to other or- 
ganic electronic devices. Based on our model, and applying a 
device simulator, we are able to predict the voltage and light 
intensity dependence of the recombination reduction factor. 
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